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Abstract. The dynamics of polymer translocation through a pore has been the subject of 
recent theoretical and experimental works. We have considered theoretical estimates and 
performed computer simulations to understand the mechanism of DNA uptake into the cell 
nucleus, a phenomenon experimentally investigated by attaching a small bead to the free end 
of the double helix and pulling this bead with the help of an optical trap. The experiments show 
that the uptake is monotonous and slows down when the remaining DNA segment becomes 
very short. Numerical and analytical studies of the entropic repulsion between the DNA 
filament and the membrane wall suggest a new interpretation of the experimental observations. 
Our results indicate that the repulsion monotonically decreases as the uptake progresses. Thus, 
the DNA is pulled in (i) either by a small force of unknown origin, and then the slowing 
down can be interpreted only statistically; (ii) or by a strong but slow ratchet mechanism, 
which would naturally explain the observed monotonicity, but then the slowing down requires 
additional explanations. Only further experiments can unambiguously distinguish between 
these two mechanisms. 
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Figure 1. Kinetic measurement of DNA uptake. (A) The nucleus being dragged after a 
3 ^tm-diameter bead, linked by a DNA molecule. The time intervals between recording the 
images are 532 s and 302 s, respectively. Note the shortening of the extended length of the 
DNA. (B) Measurements of the extended length vs. time. The data was collected from many 
measurements. Since the initial moment of association between the DNA and the nuclear pore 
was unknown, all data was fitted to a common master curve. (After M.) 



1. Introduction 



In eukaryotic cells molecular transport between the nucleus and the cytoplasm takes place 
through the nuclear pore complexes (NPCs) |[] . Ions and small proteins can diffuse through 
these large channels, but large proteins and protein-RNA complexes are transported actively 
via a signal-mediated process. DNA, on the other hand, is rarely transported under normal 
circumstances. However, this latter process might be relevant during viral infection or genetic 
therapy. In a recent elegant experiment the uptake of single DNA molecules into the nucleus 
was studied by the group of M. Elbaum [ph. In the present paper we provide a computational 
and theoretical study to understand and interpret the results of this experiment. 



2. Description of the experiment 

In the experiment double stranded bacteriophage A DNA was transported into cell-free nuclei 
of Xenopus laevis eggs. The DNA molecules were linked to polystyrene beads of 1 or 3 /im. 
Beads lying near a nucleus were trapped by optical tweezers and pulled away. In some 
cases dragging the bead pulled the whole nucleus behind it at a distance (Fig. |l]A). After 
the bead was released, it started a biased Brownian motion due to thermal fluctuations and 
its attachment to the DNA segment. Such stretches were repeated on the same bead at two 
minutes intervals, and the maximum extension decreased with each measurement (Fig. [TJb). 
This suggests that some mechanism inside the nucleus pulled the DNA inward. The remaining 
extended length vs. time curve has two main features: (i) it starts linearly, and (ii) the length 
does not go to zero, instead, it converges to about 1 fim. The uptake itself is very slow 
(compared to the time scale of biochemical processes in living organisms), it takes several 
hundred seconds for a 10 /im-long DNA. 
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3. Theoretical estimates for the entropic forces 

To explain the convergence of the remaining length vs. time curve and to give some 
estimations for the uptake force and the effective friction coefficient of the pore, Salman et. 
al ^ proposed a model in which the shrinking remaining segment of the DNA applies an 
increasing resistance against the uptake. They argue that because the entropic spring constant 
of a chain of length L is proportional to L^ 1 and the average separation of the two ends of 
the remaining DNA segment (caused by volume exclusion by the bead and the nucleus) is 
proportional to L 1 / 2 , the resulting entropic repulsion diverges as their product, L -1 / 2 . 

However, a detailed statistical analysis (described below) shows the opposite: the 
entropic repulsion imposed by the shrinking DNA segment decreases as — L^ 1 . This result 
makes it necessary to reconsider the proposed kinetics of the uptake of DNA by the nucleus. 

It has long been known that a polymer of length L, anchored at one of its ends to a wall 
(with the other end being free), pulls on the wall by the force 

k T 

Ffree ended = -Ffee ended _ (l ~ l) ~ j— (!) 

at the anchoring point, on average [0, & B]. T denotes the absolute temperature, fee is the 
Boltzmann constant, and the length independent force term F f ^ e ended > depends on the way 
of anchoring and on other properties of the polymer. 7 w 0.69 for self-avoiding and 7 = 1/2 
for non-self-avoiding polymers. Due to force balance the rest of the polymer pushes the wall 
by the same amount of force. 

Below we show that for non-self-avoiding polymers Eq. ([[]) can be obtained in a way 
much simpler than the original derivation, using a simple lattice chain model. Then we 
develop further this approach to take into account the role of the bead. 

If the persistence length of the polymer is L p , the lattice constant of the equivalent three- 
dimensional lattice chain is a = 2L p , which is called the Kuhn length [ph. Let us suppose that 
the wall is located at the x = plane, and the anchoring point is at the origin A = (0, 0, 0), 
as depicted in Fig |2|. The constraint that the chain is not allowed to pierce through the wall is 
modeled by the condition that none of the x coordinates of the chain points can be negative: 
< Xi, where < i < N and N — L/a is the number of Kuhn segments. The partition 
function Z of this system is simply the sum of all possible chain configurations that satisfy 
the above condition. 

Let us choose an arbitrary point B in the < x half-space. According to the so- 
called reflection principle, each configuration that connects A and B but does not satisfy 
the constraint, has one and only one "mirror configuration" that connects A' = (—2a, 0, 0) 
and B. Illustrated in Fig. ||a, the mirror configuration is obtained by reflecting every chain 
point from A and to the first unsatisfied point about the x = —a plane. Denoting the number 
of unconstrained configurations that connect points P and Q by J\f(P — > Q), the partition 
function can be expressed as 

Z= {^i A -> {x N ,y N ,z N )]-J\f[A' -> (x N ,y N ,z N )]\,(2) 

xn>0,vn ,zn 

where each configuration that is counted mistakenly in the first term is canceled by its mirror 
configuration in the second term. The partition function can further be written as 

Z= ^ y^i-A -» (x N ,y N ,z N )] - Af[A -> (x N + 2a,y N ,z N )]j 

%n>0,vn,zn 

Y Af[A—> (xn, yN, zjf)] (3) 

!EAr£{0,a},i/jv.ZN 
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Figure 2. Lattice chains of length L = 13a, depicted in 2D. One end of the chains is 
anchored to the membrane at point A, while the other end is either (a) free, or (b) bound to a 
mobile wall. The figures indicate how to construct the mirror configuration (light gray lines) 
of a configuration not allowed by the constraints (dark gray lines). 



fa (2d)" V(0 < x N < 2a) ps (2d)" 2a p(0) , 

where d — 3 is the dimension of the space, V(Q < xjv < 2a) is the probability that 
the chain ends between the planes x — and x — 2a in the continuum limit, and 
p(x) — (2-ko 1 )~ x I' 1 exp[— x 2 /(2a 2 )] is the probability density of the x coordinate of the 
end position with a 2 = a 2 N/d. The approximations are valid for large N. 
Thus, the partition function can be expressed as 

1/2 



r-V2 



(2d) 



N 



the free energy is 

T= -ksThi(Z), 
and the force acting on the anchored end of the chain is 



free ended 



~dL 



1 dT 

adN 



ln(2d) 
2£ n 



l fc B T 
2 L ' 



(4) 



(5) 



(6) 



in agreement with Eq. ([[]). 

In the experiments, however, the distal end of the polymer is not free, but is bound to 
a polystyrene bead. Because the bead is large compared to the remaining DNA coil when 
the uptake stalls, we can approximate it with a mobile wall. The lattice chain calculation 
can be repeated even for this system. The partition function, i.e., the number of allowed 
configurations is now 

Z= ^2 {jV[A^ (x N ,y Nl z N )]-J\f[A' ^ (x N ,y N ,z N )]- 

x N >0,y N .z N 

— N [A — > (xn + 2a, j/jv, zn)] + A/" [A' — > (xn + 2a, ijn, zn)] }(7) 
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Figure 3. The force exerted by the membrane as a function of the polymer length. The symbols 
are the results of numerical simulations with b = 50 nm, averaged over ten 2-second-long runs. 
The curves are theoretical predictions (shifted upward to fit the data set). 



where the first two terms are from Eq. (|2|), the third term - based on the reflection 
principle again (see Fig ||b) - cancels all configurations that go beyond the end 
point {xntUNtZn), and the forth term compensates for those configurations that 
have been canceled twice (in both the second and the third terms). Replacing 

TV [A' -> {x N ,y N , z N )} by TV [A -> (x N + 2a, y N , z N )], and TV [A' ->■ (x N + 2a, y N , z N )] 
by TV [A — > (a; 7v + 4a, y^, Zjy)]> and making further simplifications one arrives at 

Z= ^ N[A — > {x N ,y N ,z N )] ~ 

x N £{0,a},y N ,z N 

N[A^ {x N ,y N ,z N )\ (8) 

^jv^{2a,3a},yjv,2Af 

« (2d) N 2a {p(0) - p{2a)\ « ^(2d) N ( ^ ' 



N 



leading to the entropic force 



ln(2d) 3 k B T 

-^mobile wall ~ T— — -— — , (9) 

which indicates that a bead (or mobile wall) attached to the end of the DNA renders the 
entropic resistance even smaller. This can be understood intuitively by noticing that near 
the end of the uptake process the geometric constraints drastically reduce the available 
configuration space for the remaining segment, so it becomes more and more favorable for 
the DNA to slide through to the other side of the membrane. 

Although the reflection principle can only be applied to non-self-avoiding polymers, 
we expect a qualitatively similar behavior for self-avoiding ones, too. In this latter case the 
reduction of the configuration space should also decrease the entropic resistance. 

To verify our theoretical estimations we have performed numerical simulations. As Fig. ^] 
indicates, the numerical results for the free ended filament and for the one with the mobile wall 
are in good agreement with the relevant terms — k B T/(2L) and —3k B T/(2L) of the theoretical 
predictions, respectively. The case with the polystyrene bead, just as expected, is in between 
the above two extremes, closer to the one with the mobile wall. 
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Figure 4. Discretized polymer: the beads are connected by the rigid rods (bonds), and the 
bending elasticity is taken into account by angular springs. 



4. Simulation technique 

We modeled the double stranded DNA as a three-dimensional "bead-rod" chain [j7j], and used 
molecular dynamics simulations to follow its dynamics (see Fig. Q). Note that the "beads" of 
the bead-rod chain are artificial units and have nothing to do with the polystyrene bead in the 
real experiment. Monte Carlo methods were ruled out because they usually fail to reproduce 
the correct dynamics. 

The overdamped equations of motion of the beads are 

rRj = Vf* + Ff ad + F? etr + Ff + £j , (10) 

where T is the friction coefficient and Rj is the position of the jth bead (1 < j < N + 1). 
The bond vectors can be defined as r\, = R,-+i — Rj (for 1 < j < N). The forces on the 
r.h.s. of the equations are as follows: 

• F™ nstr is the constraint force. These forces are responsible for keeping the length of the 
rods (or bonds) fixed: |rj| = b. The constraint forces can be decomposed into the bond 
tension forces: Ff nstI = F' ens - F<™\ with F' ens = TjTj and Fjf 8 = F% a " +1 = [|]. 
The tensions Tj, represent N new variables, however, we have exactly this number of 
additional equations expressing that the bond lengths are constant: Tj ■ tj = 0. 

• F^ end is the bending force. The bending resistance is represented by angular springs 

between the adjacent bonds. The bending energy is E hend = X^ta 1 wnere K (b) 
is the bond length dependentbending modulus, 0$ = arccos(ri • ri + i/b 2 ) is the included 
angle between bonds i and i + The bending force is thus F^ end = — Vj£' bend . In the 
literature usually the bending modulus of the continuous polymer, ko, is used to calculate 
the bending energy even in a discretized polymer model. However, we realized that for 
a discretized bead-rod polymer the bending modulus should depend on the bond length 
to give correct statistics. The detailed calculations will be presented elsewhere, we show 
here only the result: k(b) w l-6/6 + 6 2 (4/7r 2 -13/24) + fe 3 (l/6-2/7r 2 ) + 6 4 /96, where 
b = b/L v is the dimensionless bond length and k(b) — k(6)/kq is the dimensionless 
bending modulus. The persistence length L p of the polymer is related to the bending 
modulus through the expression L p — ^/(k^T). A feature of this bond length 
dependent bending modulus is that k(2L p ) = 0, which reflects the known fact that a 
continuous worm-like chain with persistence length L p can be well approximated by a 
freely jointed chain with bond length b = 2L p on the large scale. 

• F™ etl is the metric or pseudo force, and is needed because the motion of a bead-rod chain 
in a thermal bath is different from the limiting case of very stiff bonds [foj |ll [l2|, |l3]] . 
The latter is a better model for real biopolymers, and the metric force turns the rigid rods 
into very stiff springs. F™ etr = —k^TVj In y\g\, where g is the metric tensor of the 
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hypersurface of the constrained N + 1 beads in the originally 3(iV + 1) dimensional 
configuration space. The explicit form of this metric tensor for a linear bead-rod chain 
can be found in Refs. |§,[l3f|. 

• F™ 1 is an optional external force acting on bead j. In the application presented in this 
paper, the polystyrene bead is attached to the end of the polymer by applying an attracting 
external force on the last polymer bead, while this bead and the wall of the nucleus exert 
repelling forces on the polymer. Depending on the situation, the bond length b may be 
too large to represent the filament contour faithfully. For example, in the simulation the 
DNA molecule goes through a 10 nm wide pore embedded in a 100 nm thick membrane 
wall, while the bod length is 100 nm. The spatial resolution of the bead-rod chain in this 
case would be very low compared to the size of the pore. For this reason we introduced 
M "virtual" beads along each rod, spaced uniformly. Their sole role is that external 
forces can act on them, and they appropriately redistribute these forces to the original 
beads at both ends of the bond they are sitting on. Hence their presence improve the 
spatial resolution, as far as external forces are concerned. We used M = 10 in the 
simulations. 

• £j represents spatially and temporally uncorrected thermal noise, obeying the 
fluctuation-dissipation relation (£)&(*')) = 2k B TY5 ]k 8(t - t'). 

In our model we neglected the excluded volume effect for the DNA filament itself, 
although we accounted for the exclusion of the filament by the polystyrene bead and the 
nuclear membrane wall. Hydrodynamic interactions were also neglected. 

In our simulations the persistence length was L p — 50 nm and the friction coefficient 
of the beads was T « irr/b, where the viscosity of the water is ?/ = 10 _3 kgs _1 m _1 and 
the temperature is T — 300 K. We used the second-order Runge-Kutta (or midpoint) method 
to solve the system of first-order ordinary equations Jl4|]. The bond length b varied between 
50 nm and 100 nm. 

One can rely on the results of a dynamical simulation only if it is numerically stable. 
Therefore, it is important to be aware of the physical processes that determine the largest 
possible time step. In our case the diffusion of the polymer beads and the relaxation of the 
angular springs pose an upper limit for the time step. The typical distance a diffusive particle 
takes during time At is Aa; = \/2D At, where the diffusion constant is D = k B T /T. To 
ensure that the typical change of the bond length is smaller than a small fraction erjb of the 
bond length, the time step cannot be larger than 

where T = T/b is the friction coefficient per unit length. Here the constraint forces were 
not taken into account. However, we expect that even with the constraint forces in effect the 
numerical error is still proportional to the typical diffusive displacement of the beads in one 
time step. 

The relaxation time of the angular springs is [ p"5[ l 

f mA 

1~AS 



6k(6) 



Apart from the prefactor, this can be easily seen from dimensional analysis. Thus, the time 
step cannot be larger than At as — £as t as either, where eas <C 1. In our simulations we 
used Ed = £as — 0.05. The maximum allowable time step is then the minimum of Atp and 
At AS - 
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5. Discussion 

A DNA filament threading through a nuclear pore is pulled by two opposing entropic forces 
from the two sides of the membrane. If the media on the two sides are identical, the length 
independent force terms cancel each other. Because the length dependent contributions (being 
proportional to k%T / L) disappear quickly for large L, there must exist some active (energy 
consuming) mechanism responsible for the pulling of the DNA. Since the integral of the length 
dependent term (from L p to the total length of the DNA) is in the order of only a few fc B T, 
this force term is not even strong enough to trap the DNA in a configuration where almost the 
entire DNA is on one side of the membrane. 

5.1. Chemical composition difference 

The simplest way to implement an active pulling mechanism is to produce a difference 
between the length independent terms of the entropic forces on the two sides of the membrane. 
This can be achieved, e.g., via a difference in the ionic strengths or the pH values of the 
two media. Even the simple lattice chain model can account for this effect if the chemical 
composition of the medium affect the persistence length L p of the DNA [see Eqs. (Q) and (0)]. 

5.2. Ratchet mechanism 

Another possibility for the pulling is based on a ratchet mechanism [|l6|]. The ratchet 
mechanism as a means of active transport was first proposed by the group of G. Oster ftl% |T8| |. 
The idea is that inside the nucleus some kind of molecule is able to bind to the DNA with a 
rate constant fc on and to detach from it with a rate constant fc ot f. When such a molecule is 
bound to the DNA at some location, the DNA cannot slide back beyond this point. 

Assuming that the possible binding sites for the molecule are separated by a distance b, 
we can estimate the effective pulling force produced by this ratchet mechanism. The grand 
canonical partition function per binding site is 2 = 1 + e~^ Eb ~ IJ '''^ kB ' r ', where Ej, is the 
binding energy and p is the chemical potential for the molecule. The probability that the site is 
unoccupied is p ff = 1/2, and the probability that it is occupied is p 0D = e~^ Eh ~^' ( fe B T )/.E. 
Thus, the partition function can be written as 2 = 1 + p on /Poff = 1 + k on /k B- The 
free energies of a binding site inside and outside the nucleus are T m — —k^T In 2 and 
!F on t = — ksT In 1 = 0, respectively. Finally, an effective pulling force can be defined as 

^ratchet = -(Fm - •7 7 out)/& = In (1 + k on /k oS ) k B T/b. (11) 

Note that this is the force necessary to stop the pulling. The pulling speed is proportional to 
this effective force only if the equilibration time of the binding of the molecules (fc on + fc ff) _1 
is shorter than the diffusion time of the DNA over the distance b. Otherwise the pulling speed 
is smaller. 

5.3. Friction in the pore 

The friction of the DNA inside the nuclear pores can significantly slow down the uptake 
process. Let us just estimate how long it would take for an entire DNA to diffuse through 
the pore if there was no friction inside. The diffusion can be envisioned as the subsequent 
diffusion of N — Ldna/o Kuhn segments (where a — 2L p is the Kuhn length). The 
diffusion time of a Kuhn segment can be approximated as tk = {l/2)a 2 ('Kr}a)/(k B T), thus 
the diffusion time of the entire DNA is about t = A 2 tk = L^ A L p (in]) / (k^T). As expected 
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Figure 5. The uptake of a DNA filament. Its length is 8 fj,m, the bond length is b = 100 nm, 
and the radius of the polystyrene bead is = 500 nm. The black dots represent occupied 
binding sites. The horizontal bar in the bottom left corner corresponds to 1 fim. The gray scale 
indicates the position in the third dimension (depth). 



for a polymer without self-exclusion, the diffusion time is proportional to the square of its 
length [[19]]. 

The diffusion time r for a 15 /im long DNA is in the order of 10 sec. Because the uptake 
occurs much slower, in the order of 1000 sec, the friction in the pore must be at least 100-fold 
greater than the hydrodynamic friction: r poie > 100L p 7T77 « 10~ 8 kg/s. Such a strong friction 
can be explained by weak binding of the DNA to the nuclear pore [p0|. 

5.4. Uptake kinetics I: weak pulling 

The experiments show that the uptake slows down in the end and the remaining length of 
the DNA converges to Lconv » 1 pm. A constant pulling force F (produced by either a 
chemical composition difference or a ratchet mechanism) corresponds to a linear potential 
U (L) = FL with a reflecting boundary at L = 0. In such a potential the expectation value 
of L at equilibrium is given by k^T jF . Thus, the measured equilibrium distance Lconv is 
compatible with an F = h&F '/Lconv w 4fN pulling force. From the average pulling speed, 
v w 2 x 10~ 8 m/s, one can also calculate the friction coefficient: r pore =f/»«2x 10~ 7 kg/s, 
which satisfies the condition given in the preceding paragraph. 

The pulling force and the friction coefficient would be even smaller if the length 
dependent term of the entropic force was also taken into account in the calculation of the 
equilibrium distance. 

We studied the uptake process by the ratchet mechanism using numerical simulations. 
We assumed that the molecules could bind only to the beads of the bead-rod chain. The 
simulations started from a configuration in which one end of the DNA had already got through 
the pore, and the first bead was not allowed to slip back through the pore. To mimic the 
experiment, a large bead of radius i?b was attached to the other, outer end of the DNA. The 
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Figure 6. The remaining length of the DNA vs. time curve. The parameters are k a = 
0.1s -1 , fc off = 0.25 s -1 , b = 100 nm, R b = 0.2 /im, and the total length of the DNA is 
4.8 fim. The curve is an average over 20 runs. 

diameter of the pore was 10 nm, the width of the wall was 100 nm, and the friction at the pore 
was neglected. Fig. || shows a sequence of such a simulation. 

Identifying the extended length of the remaining DNA by its arc length, we plotted the 
extended length vs. time curve, averaged over 20 runs. It resembles the experimental curve: 
it starts linearly and converges to a constant value, about 1 /im, which coincides with the 
experimental one. The time scale of the simulation is shorter than that of the experiments, 
because we omitted the friction inside the pore, in order to save computational time. 

The negative side of this interpretation is that the uptake is stochastic. Individual runs 
show back-slips at both the linear and the horizontal part of the curve. The stationary state of 
the system is not an energy minimum but a dynamic equilibrium between the fluctuations of 
the DNA filament and the weak pulling of the transport mechanism. Since such fluctuations 
are not apparent in the experiments, we examine another possibility for the kinetics of the 
uptake. 

5.5. Uptake kinetics II: strong but slow ratcheting 

A monotonous uptake can be achieved either by a strong continuous pulling in conjunction 
with a strong friction or by a strong but slow ratchet mechanism. In such a ratchet mechanism 
the speed of the uptake would be set by the small binding rate fc on , however, the much smaller 
unbinding rate k s <C k on would ensure that unbinding is negligibly rare and, thus, the 
ratcheting is practically irreversible. Note that the large k on /k Q ff ratio corresponds to a large 
effective pulling force. 

Due to the negligibly rare unbinding events the ratchet mechanism is also supported by 
the fact that the uptake can be stopped but never reversed in the experiments [[|] by applying 
pN forces on the polystyrene bead. A continuous pulling would be incompatible with this 
observation. 

The only phenomenon that cannot be explained by such a ratchet mechanism is the 
convergence of the remaining length to a non-zero value. There are, however, alternative 
explanations for this: 
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• the non-zero distance might be an artifact of the image processing; 

• there might be unknown structures on the nuclear envelope that prevent the polystyrene 
bead from getting in contact with the nucleus; 

• the pN forces used to extend the DNA might be sufficiently strong to deform the 
membrane by hundreds of nm-s at the vicinity of the pore [ElL 53]; 

or a combination of these. 



6. Conclusions 

Our numerical and analytical studies of the uptake of DNA filaments into the nucleus suggest 
two possible mechanisms for the uptake process. None of them can be ruled out on the basis of 
existing experimental data. One of the mechanisms explains the convergence of the remaining 
length of the DNA to a finite value, but has an intrinsically stochastic nature (back-slips occur). 
This mechanism alone cannot account for the stopping of the uptake under retarding forces in 
the pN range. The other mechanism seems more plausible, it predicts a monotonous uptake, 
but requires additional elements to explain the convergence. Only further experiments can 
unambiguously distinguish between these mechanisms, focusing on their differences. 
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